Quantum Phase Fluctuations Responsible for Pseudogap 
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The effect of ordering field phase fluctuations on the normal and superconducting properties of 
a simple 2D model with a local four-fermion attraction is studied. Neglecting the coupling be- 
tween the spin and charge degrees of freedom an analytical expression has been obtained for the 
fermion spectral function as a single integral over a simple function. From this we show that, as the 
temperature increases through the 2D critical temperature and a nontrivial damping for a phase 
correlator develops, quantum fluctuations flll the gap in the quasiparticle spectrum. Simultaneously 
the quasiparticle peaks broaden significantly above the critical temperature, resembling the observed 
pseudogap behavior in high- Tc superconductors. 
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The pseudogap, or depletion of the single particle spec- 
tral weight around the Fermi level, is the most striking 
demonstration that cuprate superconductors are not de- 
scribed by the BCS scenario of superconductivity. The 
pseudogap opens in the normal state as the temperature 
is lowered below a crossover temperature T* and extends 
over a wide range of temperatures in the underdoped 
cuprates. ARPES and the scanning tunneling spec- 
troscopy (STS) (see Refs. in provide particularly 
complete information about the pseudogap behavior and 
show a smooth crossover from the pseudo- to supercon- 
ducting (SC) gap. The transition from SC to normal be- 
havior appears to be driven by phase fluctuations Q and 
is well described by the Berezinskii-Kosterlitz-Thouless 
(BKT) theory of vortex-pair unbinding. 

There are currently many possible explanations for the 
unusual behavior of HTSC. One of these is based on the 
nearly antiferromagnetic Fermi liquid model . Another 
explanation, proposed by Anderson, relies on the separa- 
tion of the spin and charge degrees of freedom. ||^ The 
third approach, which we follow in this paper, relates the 
observed anomalies to precursor SC fluctuations, and in 
particular fluctuations in the phase of the complex order- 
ing field, as originally suggested by Emery and Kivelson 

0- 

The phase diagram for a simple microscopic 2D model 
(see (^ below) which formalizes the scenario of 
has been studied in. |^ The Green's function (GF) 
for this model was derived in ||^ using the correlator 
(exp(i6'(a;)/2) exp(— i0(O)/2)) for the phase fluctuations 
in the classical (static) approximation and neglecting the 
coupling between the spin and charge degrees of freedom. 

The associated spectral function (SF) A{ijj,\C) = 
— (l/7r)ImG(a; -I- iO, k) has also been derived analytically 
in M. Being proportional to the intensity of ARPES 
1^, the SF encodes information about the quasiparticles 
and pseudogap. The result of |^ showed that, while 
the temperature behavior of the quasiparticle peaks is 



in correspondence with experiment the gap in the 
spectrum remains unfilled. px| ] In an earlier paper pT| ] 
filling of the gap has been achieved as a result of a 
Doppler shift in the quasiparticle excitation spectrum. 
This Doppler shift originated from the semiclassical cou- 
pling of the mean field d-wave quasiparticles to the su- 
percurrents induced by classically fluctuating unbound 
vortex-antivortex pairs. It was shown in |12|] that the 
shift can be identified with the coupling of the spin and 
charge degrees of freedom. 

The purpose of the present paper is to show analyt- 
ically that filling of the gap can also result from the 
quantum (dynamical) phase fluctuations, even when the 
coupling between charge and spin degrees of freedom is 
entirely neglected. As pointed out in [ p^ the mechanism 
of the gap filling has to be understood yet even for the 
simplest attractive Hubbard model. Thus we contribute 
into the discussion ||l3| as to which mechanism or mech- 
anisms lead to pseudogap filling but within the scenario 
based on phase fluctuations, by studying the mechanism 
proposed in . The advantage of our calculation is that 
the SF is obtained as a single integral with an analyti- 
cal integrand, no numerical analytical continuation was 
necessary resulting in far greater accuracy. 

Let us consider the continuum version of the 2D attrac- 
tive Hubbard model defined by the Hamiltonian density: 

H = -^i{x) +m) M^) - Vij\{x)i^\{x)i^^{x)ij^{x), 



(1) 



where x — r,T denotes the space and imaginary time 
variables, ip^ix) is a fermion field with spin a m 
is the effective fermion mass, /x is the chemical potential, 
and V is an effective local attraction constant; we take 
h — ks — I- Clearly the Hamiltonian (|l|) is too simple 
to be adequate for systems as complex as cuprate HTSC. 
However, it has proved itself as a very convenient model 
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for both numerical, in particular Monte Carlo sim- 
ulations, and analytical approaches ||J|,|lJ,|l^ which does 
exhibit gap-like behavior above Tc (see also Refs. in the 
review |18|1 ). Moreover one may use the model to obtain 
a fully analytic treatment of the pseudogap properties, 
and apply such results to obtain a better understanding 
of more complex and less tractable models. 

The calculation of the GF is performed in Nambu vari- 
ables and the fermions are treated as composite objects, 
comprising both spin and charge parts: 

vl/t(x) - ( ^\ix) ^^{x) ) = Tt(x)cxp[-*r3^(a:)/2], (2) 

where T^(.t) is the Nambu spinor of neutral fermions. 
Substituting (||) into the standard definition of the GF 
G(x) = (vl'(a;)^''i'(0)) one obtains the GF of the charged 
(observed) fermions 

a',/3' 

as the product of the GF for the neutral fermions 
Gaf3{x) = (Tq,(x)tJj(0)) and the phase correlator 
(exp(iT36'(a;)/2) exp(— ir36'(0)/2)). For the frequency- 
momentum representation of Eq. (^ one has 

m=-Qo'' ^ ' a,f3=± 
X Dap (iuJn - il^m , k - p) , (4) 

where P± = ^(^ ± T3) are the projectors; / and 
T3 are unit and Pauli matrices; I?c(/3(«r2,„, q) is the 
Fourier transformation of the phase correlator Daf3{x) = 
{ex-p{iae{x)/2)exp{-ipe{0)/2)); ujn = {2n + l)nT and 
fin = 2m:T are respectively odd and even Matsubara 
frequencies. The GF of neutral fermions is taken in the 
mean- field approximation (see 1^,^) 



iuJni + T3^(p) - Tip 



e(p) 



P 

2m 



M (5) 



with p being a 2D vector and p = {p{x)), where p{x) 
is the modulus of a complex ordering field <^(x) = 
p(x)exp[i6'(a;)]. Note that {^{x)) = for T 7^ as 
it should be in 2D in accordance with the Colcman- 
Mermin-Wagncr-Hohcnberg theorem, while the value of 
p is allowed to be nonzero. Note also that the GF (^) 
does not contain the symmetry violating term ^ ti which 
could originate from the terms P±Q{iLOrn^)P^ since the 

correlators — = 

The representations (H), (^) for the fermion GF with 
decoupled spin and charged degrees of freedom are appro- 
priate when one can neglect the fluctuations of the mod- 
ulus p{x) and when the energy of the phase distortions 
is smaller than the energy gain due to nontrivial p. For 
the present s-wave model this means that the condition 



p^T should be satisfied, so that the modulus-phase rep- 
resentation (H) should be useful even for T close to Tbkt 
and allows one to study the evolution of the SC gap to 
the pseudogap. The region of temperatures T where the 
condition p 3> T is satisfied depends crucially on the rel- 
ative size of the pseudogap region (T* — Tbkt)/^* which 
may be reasonably large in 2D for intermediate coupling 
|19| . The pseudogap in the current work is the result of 
the low dimensionality of the system and does not need 
the existence of preformed local pairs. The nonzero value 
of p is due to average local density of Cooper pairs which 
do not have coherence above Tbkt- 

The representation (^) shows (see the explanation af- 
ter Eq. (^) that the GF for the charged fermions is de- 
fined by the correlator D{x) = Z?++(a;) = D^^(x). The 
asymptotic form of the correlator at large distances, de- 
scribing also the temporal decay of correlations, is given 
by the form 



i?«(i,r) = exp(-7i)(r/ro) 



-T/B-nJ 



exp(-r/C+(T)). (6) 



Here t is the real time, 7 is a decay constant, rg = 
(2/T)( J/if )^/^ is the scale for the algebraic decay of cor- 
relations in the SC BKT phase [T < Tbkt, where Tbkt 
is the temperature of the BKT transition), £,+ (T) is the 
phase coherence length for T > Tbkt {i+{T — > T^k^t) ~^ 
00). The constants J and K are the bare (mean-field) su- 
perfluid stiffness and compressibility, respectively which 
have been calculated in Q . 

Previously, using the representation (||), only the clas- 
sical (static f2„ = 7 = 0) fiuctuations have been consid- 
ered analytically ||]. The Fourier transform of (^ for 
this case is 



i?(ia.,q) ^ <^„,oC[q' + (i/e+)']-"/T, 



(7) 



where C = 47r(r(a)/r(l - a))(2/ro)2"-2 and a = 1 - 
T/IGttJ. For T - Tbkt the value of a ~ 1 - T/32Tbkt- 
The presence of a 7^ 1 in (^ is related to the preex- 
ponent factor (r/ro)^'^^^'"'^ inm^). This factor was not 
included in the analysis of ||Tl|,|l2|, but the treatment of 
includes this factor. 

One can now extend the analysis to the case of quan- 
tum (dynamical) phase fluctuations. We propose the fol- 
lowing generalization of (|^): 



D{ifln,q) 



Civ 



2\a 



(8) 



where v is the velocity of the Bogolyubov excitations. 
Recall that in 2D v = vfI^/2, where is the Fermi 
velocity. The asymptotic form of the retarded GF corre- 
sponding to the GF (P) is 



< 7^ - 

t^+oo. (9) 



,,2e-2 

^a-lg-*(7-\/7^-f^5+^-f^9^) ^2„2 ^ „,2 „,2c-2 
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Eq. (g) can be regarded as a convenient (for analytical 
studies) generalization of the phenomenological depen- 
dence which for nonzero 7 includes the decay of 
the phase correlations due to the presence of free vortices 
above Tbkt- One can see that for v^q^ < 7^ — w^f^^ the 
decay rate is less than 7 and it is minimal for q = Q. This 
means that, for large distances, phase fluctuations do not 
feel the pair vortices which have smaller size. 

In the present work both 7 and ^+ are phenomeno- 
logical parameters which can be derived from the theory 
of the BKT transition. Since in the SC BKT phase the 
vortices are confined, one can state that there is a 
critical slowing down of the phase fluctuations when the 
temperature approaches Tbkt, i-e. 7(T — > Tbkt) = ^■ 
In fact the detailed theory of BKT transition predicts 
that 7(T) ~ ^^^(T), where z is the dynamical exponent 
and C;\T) ~ ^Q-i exp(-6/Vr-rBKT), where b IS a pos- 
itive constant. In what follows we just restrict ourselves 
by comparing two cases: 7^0 above Tbkt and 7 = 
for T < Tbkt- 

We note that there is now experimental evidence Q 
for the vortex-pair unbinding (BKT) nature of the SC 
transition in the Bi-cuprates. We stress also that, despite 
the rather simple form of the GF, it takes into account 
the presence of vortices while the self-consistent T-matrix 
approximation (see e.g. |l3|JT^ , |2l| ) cannot describe them. 

The SF associated with GF (^) can readily be ex- 
pressed in terms of the corresponding SF for the GF (||) 
and (^) as 



00 

A(w,k) = -ilmGii(cj + iO,k) = j dio' 



1 



1 + e'^'/T 
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where these SF are 

apiuj^v) = (t^ + e(p))'5(^^^ -^^(P) -P^)sgnw, 
aB{n,Ci)^--lYnD"{^l,Ci). (11) 

TT 

The ^- function in ai? allows one to perform the angular 
integration in (|l0|). Finally integrating over momentum 
the SF can be expressed as a single integral: 



A(u;,k) 



1 



Via) 



27r^Tr(l — a) ymr. 



daj'sgn(a;' -|- uj) 



1 



X I [/(C^, k, J) (u' +UJ+ V(w'+w)2-p2j 
Xe{p + ^ [U' + LOY - p^) 

^{UJ' + Ujf - p2 ^ - V(^' + ^)2 - } , 



(12) 



where we have introduced a function /(w, k, w'): 

/(a;,k,c^') = 
vrim 

1 fuo'^ 



(.x_-a-z6)-"T(l,a;l;^±-^ 
'2 a + lb — X- 



x± 



70; 



(13) 



For 7 = the expression ([l^ can be rewritten in the 
following form 



A(a;,k) = 



1 



r(a) 



Q-1 



xsgn(aj' -I- uj) 



27r2Tr(l-a) \mrl 
1 1 



dw'sgntj' 



^{U' + Lof - p- 



nia - a;_)-"T ( i, a; 1; 0{a - x+) 

+ (x+ - x^)-'/'(a - x_)i/2--s(l , 1 _ a) 



xF 
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a; ) 0{x-^ — a)0{a — x_) 



X (c^' + + ^{oj' + ojy - p^^ 6{^i + ^{uj' +uf - p'^) 
+ y{uj' + ioY - p2 ^^{uj' + u:Y - p2) } (14) 

which describes the case of the non-damped dynamical 
phase fluctuations. While Eq. ( p^ ) appears to be more 
complicated than the more general Eq. (^2|), it proves 
useful for making a general statement about the gap fill- 
ing- 

Expressions ( |l2| ) and (|l4| ) present the main result of 
the paper for the SF in the case of quantum phase fluc- 
tuations. Before proceeding to the discussion of their nu- 
merical integration, we recap briefly the case of classical 
phase fluctuations. The static SF, obtained in |^, is re- 
produced from (|l^) after changing the variable lo' vlu' 
and taking formally the limit v —^ 0. It was discussed 
in detail in Q and we would like to stress here only two 
main points. The first is that the SF is identically zero 
inside the gap {Aci{LU,k) = for \lu\ < p) The 
second is that in addition to the usual quasiparticle peaks 
it reveals extra peaks at u = ±p. 

Let us now return to the discussion of the results of 
numerical integration based on Eqs. (|l2|) and (|lj) which 
are shown in Fig. |l|. 

a) For T < Tbkt there are two highly pronounced quasi- 
particle peaks at w = ±i?(k) = ±y/^^k) + p^ . Since be- 
low Tbkt both and 7 are zero, the width of the peaks 
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is almost entirely controlled by a whose deviation from 1 
gives the non- Fermi liquid behavior of the GF ^ . Since 
the value of a is very close to 1 this non- Fermi liquid be- 
havior is hardly distinguishable from ordinary widening 
of the quasiparticle peaks due to damping. Furthermore, 
because a ^ 1 as T ^ the width of the peaks is tem- 
perature dependent so that the peaks become sharper as 
T decreases resembling the data of ARPES for the anti- 
node direction. |^ 

We note also that k ^ kp iw Fig. |l| so that the quasi- 
particle peaks are not equal to each other in contrast to 
the case of symmetrized ARPES data [|l| . The reason for 
choosing k ^ kp 'va. Fig. |^ was to prove that the extra 
peaks present in the SF calculated for the classical phase 
fluctuations are now absent. 

The gap in the SF remains almost unfilled and has 
"U"-like shape. We stress, however, that in contrast to 
the static case |^ the SF ( p^ is nonzero even for < p 
and this is evidently related to the presence of the quan- 
tum (dynamical) phase fluctuations, i.e. the terms with 
ri„ ^ 0. In particular, using Eq. ( [l4| ) one can check an- 
alytically that even for 7 = that due to these terms 
A{uj = 0,k = k^) 0. We estimated the temper- 
ature Tci of the quantum to classical crossover in the 
way similar to that of psf . For the low carrier densities 
this gives Td ~ Tbkt, while for the high carrier den- 
sities, where the pseudogap phase shrinks [|j we obtain 
Tci ^ p{T = 0) ^ T*. Thus, the crossover tempera- 
ture is always too large (in the absence of dissipation) 
for classical fluctuations to play a significant role at low 
T. 

b) For T > Tbkt the quasiparticle peaks at w « i^(k) 
are far less prominent. This is caused by the fact that 
^+ is now nonzero due to the presence of vortices. In- 
creasing the temperature further decreases the value of 
^+ so that the quasiparticle peaks become even less pro- 
nounced. This behavior of the quasiparticle peaks repro- 
duces qualitatively the ARPES studies for the anti- 
node direction which show clearly that the quasiparticle 
SF broadens dramatically when passing from the SC to 
normal state. The width of the quasiparticle peaks is 
controlled primarily by ^+ above Tbkt since this width 
remains practically constant as 7 changes. Since ^+ is the 
phase coherence length, the current model supports the 
premise that the quasiparticle peaks are related to the 
phase coherence, not to the energy gap. The deviation 
of a from 1 becomes, however, sizable at T > Tbkt- 

c) As discussed above, for T > Tbkt the value of 7 is 
nonzero and increases with increasing temperature. This 
increase of 7 (typically up to values of the order of 0.5/i), 
together with the decrease of , causes filling of the gap 
and changes its shape from "U"-like to "V"-like. In other 
words, the quasiparticle peaks grow "shoulders" which 
eventually fill the gap. We stress, however, that the in- 
crease of 7 in fact diminishes the effect of the quantum 
phase fluctuations pushing Td down ||23|| which is seen 
from the reducing of A{uj = 0,k = \ip) as 7 grows (see 
Fig. 1, where the bottom parts of the curves cross each 



other), but the abovementioned "shoulders" simultane- 
ously emerge. Thus the filling of the gap at T > Tbkt 
in the present model is caused by quantum phase fluctu- 
ations in the presence of nonzero dissipation because the 
damping 7 contributes only when Qn 7^ 0. 
d) Due to the smooth dependence of and 7 on T as 
the temperature passes through Tbkt there is no sharp 
transition at Tbkt in agreement with experiment ||^,^ . A 
similar filling of the gap was obtained in jl^ for 7 — O.bfi 
where the correlator (|^) was used for the numerical com- 
putation of the self-energy of the fermions and the subse- 
quent extraction of the SF from the fermion GF. A recent 
Monte Carlo simulation also shows similar behav- 
ior for the quasiparticle peaks and the filling of the gap. 
However it is the analytical character of the present work, 
which relies on the explicit introduction of the charge and 
spin degrees of freedom (H) for the Nambu spinors, that 
enables us to unambiguously state the correspondence 
between the parameters of the model and the observed 
features of the SF (|l|). 

In spite of some similarities between the results ob- 
served and those just obtained, the latter are only il- 
lustrative since we have considered a model with non- 
retarded s-wave attraction. However, it is likely that for 
d-wave pairing, the properties obtained can be used for 
the description of the systems in the anti-node direction 
on the Fermi surface. It is however essential to consider 
fluctuations in the modulus of the order parameter to 
extend the analysis to the nodal directions on the Fermi 
surface. The value of 7 = 0.5/i which results in a filled 
gap appears to be too large to be due to the vortex- 
vortex interaction. This leads one to the conclusion that 
the mechanism considered here for gap filling may well 
not be the only possible mechanism and that other inter- 
actions which lead to the filling of the gap above Tc, see 
]ll[-p^, are also important. 

In summary we studied the effect of the fluctuations in 
the phase of the complex ordering field on the properties 
of a 2D system with four-fermion attraction. The fermion 
SF has been given as a single integral over a function 
known in closed form. Through the use of analytical 
techniques, we have been able to demonstrate that the 
quantum phase fluctuations in the presence of dissipation 
lead to the filling of the pseudogap even if one ignores the 
spin-charge coupling. 
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FIG. 1. The spectral function A{uj, k) as a function of uj in 
units of the zero temperature SC gap A for k > kp taken at 
two different temperatures, Ti < Tbkt < Tb. 
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